WMAP Constraints on Decaying Cold Dark Matter 
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We re-formulate cosmological perturbations in the decaying cold dark matter model, and calculate 
cosmological microwave background (CMB) anisotropics. By comparing our theoretical predictions 
with recent observational data from the Wilkinson Microwave Anisotropy Probe (WMAP), we derive 
a new bound on the abundance and lifetime of decaying dark matter particles. We show that the 
data of WMAP alone do not prefer the decaying cold dark matter model; the lifetime is constrained 
to F"^ > 123 Gyr at 68% C.L. (52 Gyr at 95.4% C.L.) when cold dark matter consists only of 
such decaying particles. We also consider a more general case in which cold dark matter consists 
of both stable and decaying particles, and show that the constraint generalizes down to noDM^i'^ < 
-0.5(F-VlGyr)-^ + 0.12 for F"^ > 5Gyr at 95.4% C.L.. These bounds are robust and widely 
applicable, because they are derived only from gravitational effects and therefore do not rely on the 
details of decay channels or decay products. 



Introduction. — The existence of cold dark matter 
(CDM) is now widely accepted from numerous kinds 
of astronomical phenomena, such as rotation curves in 
galaxies, anisotropies in cosmic microwave background 
(CMB), and X-ray emitting clusters of galaxies combined 
with the low cosmic baryon density predicted by big bang 
nucleosynthesis. The nature of CDM, however, is still one 
of the biggest mysteries in cosmology. Indeed, some dis- 
crepancies on galactic and sub-galactic scales in standard 
CDM cosmology have stimulated numerous proposals to 
modify the standard CDM model 0. 

The decaying cold dark matter (DCDM) model is one 
of such proposals to solve discrepancies in the standard 
CDM model. For instance, Cen proposed the DCDM 
model to solve both over-concentration problem of the 
dark matter halos and over-production problem of small 
dwarf galaxies The authors also showed in previ- 
ous papers '3| that introducing DCDM can improve the 
fits of observational data sets of Type la SN, mass-to- 
light ratios and X-ray gas fraction of clusters of galaxies, 
and evolution of cluster abundance. The DCDM model 
may be also useful to solve other important problems 
such as ultra-high energy cosmic rays above the Greisen- 
Zatsepin-Kuzmin cut-off of the spectrum Q and the early 
reionization of the universe f^. 

Although theoretical candidates for DCDM have been 
proposed, their predictions for lifetime of decaying par- 
ticles cover a large range of values, from 10~^ Gyr to 
10^^ Gyr 1^. This indicates that it is quite important to 
constrain (possible) lifetime of dark matter particles from 
astronomical observations. One popular way to constrain 
the lifetime is to use decay products, e.g., diffuse gamma 
ray background observations 0- However, realistic sim- 
ulation which takes all energy dissipation processes into 
account showed that even the particles with lifetime as 
short as a few times of the age of the universe still are not 
ruled out by recent observations Moreover, the re- 



sults are dependent on the assumed decay channel while 
we do not know what dark matter is. 

In this paper, we derive cosmological constraints on 
the abundance and lifetime of decaying dark matter par- 
ticles from CMB anisotropies. We re-formulate cosmo- 
logical perturbations in the decaying cold dark matter 
model and compare them with accurate measurements 
by the Wilkinson Microwave Anisotropy Probe (WMAP) 
jSl . Here we assume only that dark matter particles decay 
into relativistic particles, whatever they may be. Thus 
our constraints are different from those derived from de- 
cay products such as diffuse gamma ray background ob- 
servations which need to assume the details of the decay 
products or decay channels. Throughout this paper we 
concentrate our attention on a flat universe. 

DCDM cosmology: Basic Equations. — To make the- 
oretical predictions of CMB anisotropies in the DCDM 
model, first we must formulate cosmological perturbation 
theory. Similar works can be found in references [1(1 
in which they discussed decaying neutrinos. The present 
work differs from them in that we consider decays of cold 
dark matter, and we treat both the abundance and the 
lifetime of decaying particle as free parameters. 

The equations of background energy densities for 
DCDM particles (DDM) and its relativistic daughter par- 
ticles (DR) are given by 



PDDM — — 3-ffpDDM — arpDDM 
PDR = -4-ff/ODR + arpDDM ■ 
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Here dot denotes conformal time derivative and F is the 
decay width of the dark matter. The equation of state 
parameters w = P/ p are ioddm = , and wdr = 1/3 , 
respectively. We also define effective equation of state. 
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which are defined by the evolution of energy density, p 
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-3H{1+Wes)pi. 

Now let us turn to the cosmological perturbations. In 
the conformal Newtonian gauge, the hue element is given 
by 

ds^ = a(T)2 [-(1 + 2ijj)dT^ + (1 + 2<i>)S,jdx'dx^] . (4) 

where a(r) is cosmic scale factor, ^ and $ are perturba- 
tions around flat Friedmann-Robertson- Walker metric. 

To derive the equation of DCDM we start from the 
Boltzmann equation for non-relativistic matter with de- 
caying term and take its zeroth and first moments by 
momentum integration: 

n + di{nv') + 3{H + = -aVn , (5) 
drijiv-') + AHnv^ + ndjtjj — —aVnv-' . (6) 

We decompose n(r, x*) = no(T)(l + 6{t,x^)) and use 
gauge invariant density perturbations defined by = 
6 + 3(1 + Wcs)^- Then we have the following equations 
for DCDM: 
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here we have defined v = i{kj/k)v^ . 

Once we have equations for DCDM we have those 
for its daughter radiation (DR) from the consequence of 
energy-momentum conservation: 

+ 3w^^<i> + 3(1 - _ 3^D/ )H<^> 
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3{l + w^^^^)<f) , (9) 



VDR - -l)HvuR + kib-^kIloR + ^Dl^^ 



-(1 + Wcff + aVvuBM —. 

4 V 4/ODR 



(10) 



Further, we need Boltzmann hierarchy for I > 2 mo- 
ments of daughter radiation, where / stands for multipole 
moment in a Legendre expansion of perturbed distribu- 
tion function. Again we begin with the Boltzmann equa- 
tion for relativistic particle |0| , 
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where q = qh is the comoving 3-momentum with 
n^rii — 1. We wrote the distribution function of daugh- 
ter radiation with background distribution and pertur- 
bation around it as /^^ (a;*, g, , r) = /^^(g,r)(l -t- 
'^^^{x'^ ,q,nj,T)). We should note that, unlike the stan- 
dard CDM, f^^ is time dependent in the DCDM model. 



To describe the decay process, let us consider the collision 
term: 



dl 
dr 
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where /°dm = /DDM(^ ^ ^ddM) ^^^^ ^-^e 
distribution function and mass of DCDM particles, re- 
spectively. Then we obtain 



-I- ikpF^^ -f 4($ + iknij}) 
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where /i = k ■ fi, and we defined _F^^ (fc,n, r) = 
/ q^dq (?/o°^*°^// q^dq qf^^ . The two terms in r.h.s 
of Eq. ifT^ arc new ones in the DCDM model. The first 
term comes from the fact that evolution of background 
energy density of daughter radiation is different from de- 
coupled radiations which diminish as a~'^. The second 
term corresponds to the flow from the dark matter to 
daughter radiation in first order perturbation. This is 
complicated to compute in general but it can be 
easily described in a fluid approximation (see Eq. @ 
and Ea. pU|) for I ~ and 1 = 1, respectively). For 
higher multipoles (l > 2), this term vanishes since the 
perturbation in the dark matter (^i'^'^^) does have only 
I = and 1 perturbations which correspond to density 
and velocity perturbations, respectively. Finally, we have 
hierarchy for daughter radiation (/ > 2), 
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where Mi is the coefficient in a Legendre expansion, 
If™ + m = M = J2i-iyMi{T,k)Pi{p) . These co- 
efficients are related with above perturbation variables 
by D'^^ = 4Mo, VBR = Ml, and Hdr = (12/5)M2. 

CMB constraints and Discussions. — We are now ready 
to compute theoretical CMB spectrum in DCDM cos- 
mology. We first consider a simple DCDM cosmology in 
which CDM component consists only of decaying parti- 
cles, and constrain the abundance and lifetime of DCDM 
particles from observed CMB anisotropics of the WMAP. 
Next, as a generalization of DCDM cosmology, we inves- 
tigate the universe in which CDM consists both of stable 
(standard) and decaying particles. 

The likelihood functions we calculate are given by 
Verde et al.0|. To include the decay of dark matter 
particles described above, we calculate theoretical angu- 
lar power spectrum of CMB fluctuation {l{l + 1]Ci/{2'k)) 
using a modified Boltzmann code of CAMB which 
is based on a line of sight integration approach |l5l | . 

An illustration of CMB power spectrum in the DCDM 
models is shown in Fig. ^ The decay of CDM parti- 
cles modify the CMB spectrum mainly in three ways. 
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FIG. 1: CMB angular power spectrum with and without 
decay of CDM particles, and observational data from the 
WMAP. The solid line corresponds to the best-fit standard 
CDM model. The other hues are DCDM models with T'^ = 
13 (dash), 23 (dotted), and 43 Gyr { dash- dotted) . For each 
lines, all cosmological parameters (except for lifetime) are 
fixed to WMAP optimal values to demonstrate how decays 
of dark matter particles modify the angular power spectrum. 

First, the modified evolution of dark matter changes the 
expansion history of the universe and generally causes 
the shorter look back time to the photon last scattering 
surfaces. This pushes the acoustic peaks in the CMB 
spectrum toward the larger scales (smaller Second, 
decays of dark matter particles mean that there is more 
dark matter at earlier times. This results in smaller 
anisotropics around the first acoustic peak. Third, decays 
of gravitational potential at later epochs due to decays 
of dark matter particles lead to a larger late integrated 
Sachs- Wolfe (LISW) effect. The LISW significantly en- 
hances CMB anisotropics at low multipoles. Thus, the 
amplitudes and locations of the peaks in the power spec- 
trum of microwave background fluctuations ^l(f\ can in 
principle be used to constrain the DCDM model. 

It is well known, however, that there are other cosmo- 
logical parameters which also modify the CMB spectrum. 
Therefore, we have to generate the full probability dis- 
tribution function and marginalize over nuisance param- 
eters to obtain the constraint on parameter(s) which we 
are interested in. To realize this, we followed the Markov 
Chain Monte Carlo approach [13| and explore the likeli- 
hood in seven dimensional parameter space i.e., five stan- 
dard parameters, flhh^ (baryon density), h (Hubble pa- 
rameter) , Zj-e (reionization redshift) , Ug (power spectrum 
index). As (overall amplitude), and two DCDM parame- 
ters, J7ddm/i^ (DCDM density) and F (decay width). 

Our results are shown in Figs. |2]and|21 An interesting 
point is that the parameter of the DCDM model, F, dose 
not degenerate with the other cosmological parameters 
very much as one can see in Fig. |21 This means that the 
change in the CMB spectrum from F cannot be mimicked 
by other standard parameters. This is the reason why 
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FIG. 2: Contours of constant relative probabilities in two 
dimensional parameter planes. Lines correspond to 68% and 
95.4% confidence limits. 
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FIG. 3: Marginalized probability distribution of the decay 
parameter F in Gyr~^. Confidence limits of 68% and 95.4% 
are shown by shaded regions. 



CMB can put strong constraint on the lifetime of DCDM 
particles. Figure 13 shows marginalized likelihood of the 
lifetime of DCDM particles: the constraint is F^^ > 52 
Gyr at 95.4% C.L., and T'^ > 123 Gyr at 68% C.L.. 
We do not find any signal to prefer the DCDM model to 
standard CDM cosmology. 

To obtain more general constraints on the DCDM 
models, next we consider a cosmological model in which 
CDM consists of two types of particles, i.e., one is stable 
and the other is unstable. We performed the same CMB 
likelihood analysis described above with another cosmo- 
logical parameter, J7nDM, which stands for the current 
density of stable dark matter in critical density units. 
Our results are shown in Fig. 01 The results are easily 
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FIG. 4: The same as Fig|21 but for cosmology in which dark 
matter consists of two species (stable and decaying dark mat- 
ters). The densities of stable and decaying dark matters are 
denoted by f2nDM and JIddm, respectively. 

understood: first, the universe should have either CDM 
or DCDM component f^cDM^i^ « 0.1, though CDM are 
slightly favored compared with DCDM, (left of Fig. 0J; 
second, if the dark matter mainly consists of unstable 
particle, its lifetime is tightly constrained (right of Fig. 
21) as shown in Figs. |21and|31 However, as the fraction 
of decaying components becomes smaller, constraints on 
F"^ becomes weaker. We find that the lifetime of unsta- 
ble CDM particle is roughly constrained to 

^BDMh^ < -0.5 (^) ^ + 0.12 for T'^ > 5Gyr 
^DDuh^ < -0.03 (th^) ^ + 0.026 for F-^ < 5Gyr 

at 95.4% C.L.. Bounds on DCDM would not be improved 
so much by future experiments since DCDM affects CMB 
fluctuation only in rather large scales (see Fig. QJ, and er- 
rors in these scales are already limited by cosmic variance. 
In this sense our results are robust and the tightest con- 
straints on DCDM cosmology from CMB anisotropics. 

Before conclusion, we briefly mention the reionization 
history in DCDM models. If the energy of decaying par- 
ticles is transformed into high energy photons, they may 
become a source of ionizing photons and lead an early 
reionization of the universe. This results in the suppres- 
sion of acoustic peaks in the CMB spectrum 18], which 
is similar to the purely gravitational effect considered 
here. Therefore, constraints on the decay rate might be 
stronger when we include the ionizing photons produced 
from the decays. However, this process depends on the 
details of the decay process, such as the power index of 
the injected spectrum and decay modes to a certain ex- 
tent. On the other hand, our constraints are based only 
on the minimal assumption that dark matter particles 
decay into relativistic particles. Thus our result is the 
most conservative one. 

Conclusion.- We showed that even the current CMB 
data alone put strong constraint on the lifetime of cold 
dark matter to T'^ > 123 Gyr at 68% C.L. (52 Gyr 
at 95.4%). By considering both stable and decay- 



ing dark matters, the constraint on DCDM with life- 
time larger than 5 Gyr generalizes down to Oddm^^ ^ 
-0.5(F- VlGyr)-! + 0.12 at 95.4% C.L.. The results are 
based on the minimal assumption, free from details of 
decay products and their propagation, and thus widely 
applicable. 
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